proc mixed data=WORK.CHINA covtest CL;
    class YEAR COHORT5;
    model CHINA = AGE5 / solution CL;
    random YEAR COHORT5 / solution; 
    title "Table 1 Model 1";
run;

proc mixed data=WORK.CHINA covtest CL;  
    class YEAR COHORT5;
    model CHINA = AGE5 BIEDU REGION_dummy1 REGION_dummy2 REGION_dummy3 REGION_dummy4 REGION_dummy5 REGION_dummy6 SEX_dummy2  / solution CL;
    random YEAR COHORT5 / solution;
    title "Table 1 Model 2";
run;

proc mixed data=WORK.CHINA covtest CL;  
    class YEAR COHORT10;
    model CHINA = AGE5 / solution CL;
    random YEAR COHORT10 / solution;
    title "Appendix Table A2 Model 1";
run;

proc mixed data=WORK.CHINA covtest CL; 
    class YEAR COHORT10;
    model CHINA = AGE5 BIEDU REGION_dummy1 REGION_dummy2 REGION_dummy3 REGION_dummy4 REGION_dummy5 REGION_dummy6 SEX_dummy2/ solution CL;
    random YEAR COHORT10 / solution;
    title "Appendix Table A2 Model 2";
run;

proc mixed data=WORK.RUSSIA covtest CL;
    class YEAR COHORT5;
    model RUSSIA = AGE5 / solution CL;
    random YEAR COHORT5 / solution;
    title "Table 1 Model 3";

proc mixed data=WORK.RUSSIA covtest CL;
    class YEAR COHORT5;
    model RUSSIA = AGE5 BIEDU REGION_dummy1 REGION_dummy2 REGION_dummy3 REGION_dummy4 REGION_dummy5 REGION_dummy6 SEX_dummy2  / solution CL;
    random YEAR COHORT5 / solution;
    title "Table 1 Model 4";
run;

proc mixed data=WORK.RUSSIA covtest CL;
    class YEAR COHORT10;
    model RUSSIA = AGE5 / solution CL;
    random YEAR COHORT10 / solution;
    title "Appendix Table A2 Model 3";

proc mixed data=WORK.RUSSIA covtest CL;
    class YEAR COHORT10;
    model RUSSIA = AGE5 BIEDU REGION_dummy1 REGION_dummy2 REGION_dummy3 REGION_dummy4 REGION_dummy5 REGION_dummy6 SEX_dummy2  / solution CL;
    random YEAR COHORT10 / solution;
    title "Appendix Table A2 Model 4";
run;

proc mixed data=WORK.US covtest CL;
    where YEAR > 2006;  
    class YEAR COHORT5;
    model US = AGE5 / solution CL;
    random YEAR COHORT5 / solution;
    title "Table 1 Model 5";
run;


proc mixed data=WORK.US covtest CL;
    where YEAR > 2006;  
    class YEAR COHORT5;
    model US = AGE5 BIEDU REGION_dummy1 REGION_dummy2 REGION_dummy3 REGION_dummy4 REGION_dummy5 REGION_dummy6 SEX_dummy2  / solution CL;
    random YEAR COHORT5 / solution;
    title "Table 1 Model 6";
run;

proc mixed data=WORK.US covtest CL;
    where YEAR > 2006;  
    class YEAR COHORT10;
    model US = AGE5 / solution CL;
    random YEAR COHORT10 / solution;
    title "Appendix Table A2 Model 5";
run;

proc mixed data=WORK.US covtest CL;
    where YEAR > 2006;  
    class YEAR COHORT10;
    model US = AGE5 BIEDU REGION_dummy1 REGION_dummy2 REGION_dummy3 REGION_dummy4 REGION_dummy5 REGION_dummy6 SEX_dummy2/ solution CL;
    random YEAR COHORT10 / solution;
    title "Appendix Table A2 Model 6";
run;
